Negative differential magneto-resistance in ferromagnetic wires with domain walls 
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A domain wall in a ferromagnetic one-dimensional nanowire experiences current induced motion 
due to its coupling with the conduction electrons. When the current is not sufficient to drive the 
domain wall through the wire, or it is confined to a perpendicular layer, it nonetheless experiences 
oscillatory motion. In turn, this oscillatory motion of the domain wall can couple resonantly with 
the electrons in the system affecting the transport properties further. We investigate the effect of 
the coupling between these domain wall modes and the current electrons on the transport properties 
of the system and show that such a system demonstrates negative differential magnetoresistance due 
to the resonant coupling with the low-lying modes of the domain wall motion. 

PACS numbers: 75.60.Ch, 75.78.-n, 85.70.-w, 75.47.-m 
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I. INTRODUCTION 

When a current is passed through a domain wall (DW) 
in a ferromagnetic wire then the spin orientation of the 
current electrons is changed as they pass through the 
DW, both due to adiabatic changes and scattering from 
the noncollinear magnetization. Conservation of angu- 
lar momentum tells us that this must be compensated. 
Some of the angular momentum is transferred to the lat- 
tice and some is compensated by a change in the an- 
gular momentum of the magnetization comprising the 
DW. This change in angular momentum effectively sets 
the wall in motion, which can be described using, for in- 
stance, the Landau-Lifschitz-Gilbcrt (LLG) equation^ - — 
The details of the carrier-DW coupling, and the subse- 
quent magnetization dynamics of the DW set into motion 
by the applied current have recently been the subject of 
a large volume of work^ which may also have possible 
applications to memory device technology^ The interac- 
tion of the current electrons with the domain wall can 
even induce interaction effects between the DWs»^r— 

In this paper we wish to address effects which occur 
when the domain wall is pinned. If the current is not large 
enough to cause lateral motion of the DW through the 
wire it is nevertheless possible for the DW to experience 
oscillatory motioni 7 -^ - — Such oscillations of the DW 
have been experimentally verifiecU^ and can be thought 
of as the zeroth mode excitations of the DW. In this 
paper we describe the effect that the coupling between 
such DW modes and the conduction electrons has on 
the current through the wire. The resonant coupling be- 
tween the conduction band and the mode of the domain 
wall's motion leads to a negative differential magnetore- 
sistance (NDMR). Negative differential resistance (NDR) 
is a well known phenomeno n 16 ' 17 that has been stud- 
ied extensively, e.g. for semiconductor-based resonant- 
tunneling diodes, for organic semiconductor spin valve 
systems , 18 ' 19 and for double quantum dot set-ups.— NDR 
has found diverse applications, e.g. for microwave signal 



amplifications,— in feedback oscillators,— in frequency 
mixers*^ and others. To our knowledge, however, a DW 
negative differential magneto-resistance has not been in- 
vestigated yet. In view of the NDR applications men- 
tioned above and the fact that domain wall logics is well- 
established by now,— it is timely to consider NDMR in 
noncollinear magnetic structures. 

The model we will consider is also appropriate for con- 
sidering a three part set up consisting of two ferromag- 
netic wires of opposite magnetic orientation with a strip 
of ferromagnetic material between them. The middle 
layer will have a magnetization perpendicular to the two 
wires and as a current is passed through a spin torque 
also acts on this thin layer causing its magnetic orien- 
tation to move^ Such structures are also experimentally 
promising from our point of view. 

In addition to the coupling between the DW mode 
and the electrons there arc contributions to the resis- 
tance from scattering from the domain wal k 25-28 These 
contributions to the current are included exactly within 
our model and are naturally the cause of the DW's mo- 
tion. The scattering of the electrons from the DW leads 
to a charge and spin build up in its vicinity and Friedel 
oscillations in the spin-dependent density of the carrier 
electrons. This build up of spin enhances the coupling 
between the DW motion and the electrons. 

Our method is to first solve the effect of the DW and 
its motion on the charge build up in the system. The 
motion of the domain wall is modeled by the DW oscil- 
latory mode coupled to the electrons in the vicinity of 
the DW. 11 ' 12 ' 18 ' 19 This oscillatory mode is described by 
a free energy for the classical magnetization. 



II. MODEL 

We start from a model for the classical inhomogeneous 
magnetization that we describe by a time dependent unit 
vector field n(r , t) since we assume the longitudinal dy- 
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namics are energetically forbidden. We include the ef- 
fects of the anisotropy and exchange and the coupling of 
n(r, t) with the conduction electrons with a Kondo-type 
coupling strength J. T is an applied field modeling the 
torque caused by the conduction electrons and S is the 
spin density of the conduction electrons. The magneti- 
zation free energy then rcadsi 



J" = 



Ctik 



dridr k 



T-JS- Knpf3 



(1) 



/o is the homogeneous part of the exchange energy, on k 
is a tensor describing the inhomogencous part of the ex- 
change energy and is of order ~ T c /a where T c is the 
Curie temperature and a is the lattice spacing. K gives 
the anisotropy, taken in the $ direction. 

For the conduction electrons we have, in addition to 
the coupling to this bulk magnetization, the Hamiltonian 
(we set H = 1 throughout) 



H, 



J dV^ct^-^r) 



(2) 



where e is the single particle dispersion measured with 
respect to the Fermi level [i and c a (f) are the carrier 
field operators with a spin index a. The total system is 
then described by T + H e . 

Upon perturbation, the magnetization n(r, t) under- 
goes some periodic motion in the vicinity of the pinned 
DW. The coupling between the motion of the DW and the 
conduction electrons can be treated as a bosonic mode 
of the system. Thus, integrating over the spatial direc- 
tions perpendicular to the wire, we can write down the 
following one dimensional Hamiltonian: 

H = H + J dz^ct(z)g„„>(z)(<i>i+4)c<7'(z) 

+ ^eir4>i<f>i, (3) 

i 

where ei r are the energies of the DW's resonant motion, 
9aa'(z) = g<T<r>S(z) = ga^ a ,5{z) is the coupling between 
the DW mode and the electrons, and <j> models the DW 
mode. Ho is the "uncoupled" part of the Hamiltonian. 
It is important, for the possibility of observing a negative 
differential resistance effect, that the energy of the mode 
of motion of the DW is of the order of the carrier elec- 
trons' energy. This requires cither rather narrow domain 
walls or a suitable layer of magnetic noncollincarity. Here 
we consider a system with a relatively low carrier density 
such that their Fermi wavelength Xp > L, where L is the 
DW width. Therefore it is natural to model the system 
with sharp domain walls. This has the further advan- 
tage of analytic tractability. The carrier scattering from 
a DW has been addressed previously^ From the trans- 
mission and the reflection coefficients we find the charge 
build up in the system due to the DW in the presence 



of an applied bias. In this case we deal with a standard 
s — d Hamiltonian 



H = I dzJ2cl(z)[e- fi]c a {z) 



(4) 



Where 



/oo 
dzJn x (z) w JirL and 
- oo 

lim Jn z (z) — > — Jsgn(z). 



(5) 
(0) 



Without loss of generality we have orientated our domain 
wall in the x direction. We have implicitly assumed a 
separation of time scales which allows us to treat the 
domain wall as adiabatic on the time scale relevant for 
the dynamics of the conduction electrons. 



III. CALCULATION 

The electronic Green's function is 
iG{z,z';t,t') = I 
with the action, S — + S$ + , given by 



DxpDxpD<l>D<H) z {t)%l) z ,{t'y b (7) 



/ dtd^y^ 1p<r(z,t) idtSaa' - Hq 
Jc 

- j didz^i/v(z,t) [g<Ta'{z) {<!> + 4>)] VV'OM) 
di0 [id t - e ir ] <j>. (8) 

S v ' 



For convenience c is the full interaction time contour^ 
One can fully integrate out the bosonic degree of free- 
dom, and then expand perturbativcly in the coupling g. 
Switching to the Keldysh representatio n 3 1 ' 32 and intro- 
ducing the following retarded/advanced bosonic Green's 
functions, P r / a {uj) = D R l A (u) + D r ' a {-uj), gives the 
first order correction in perturbation theory: 



G in (z, z'; e) rj G° in (z, z 1 ; e) 



(9) 



+i 



dzi dz 2 g r7l a 2 {zx)g a3rTi (z 2 hj k 7f m 



P sp (uj)G% (z,z 1 ;s)G° kl (z 2 ,z 3 ;e-io)G° m n(z 4 ,z / ;e) 
-P s P(0)G° lm (z,z 2 ;e)G° ln (z 2 , z'; e)G° kj (z u zx\ui) 



CTCT3 



We use Latin indices for the Keldysh matrix indices. 
Summation over repeated indices is implied. The emis- 
sion/absorption tensors are 7^ = and ~ff k = af k . 



3 



o20 



^2 




O 
O 




0.02 0.03 
eV/e^ 



0.05 



FIG. 1. (Color online) The x- and the z component (inset) of 
the spin density S(z), around the domain wall (of width L) 
caused by scattering of electrons from the wall. 



The first step is to calculate the spin density of the 
conduction electrons which couples to the magnetiza- 
tion, determined by Eq. (j4]). It is helpful to first de- 
compose into the scattering states i/) n (z,(r) labeled by 
n± = {±k £l a}££ (a is the spin and ±fc E the momentum 
of the incoming electron.) Then the spin density becomes 
S(z) = S L {z) + Sr(z) where 



Sl,r{z) 



E 



deis a {s)ip* ± (z, a)ip n± (z, a')a a 



(10) 



The L and R indices naturally refer to electrons incident 
from the left and right, which are held at a potential drop 
of eV symmetrically across the Fermi energy ef- Va{s) is 
the density of states of the carriers. In the region around 
the DW the transpose component of the spin density is 
enhanced, see Fig.Q] Density fluctuations are also clearly 
visible. 

The uncoupled Green's function in the scattering basis 
becomes 



G aa ,(z,z';e) = J2 



jj n {z,a)ipl(z',cr') 
£ - £ n + [J. c + iS sgn(e) 



using the shorthand 



E 

= {L,R} a 



E / de n v c a (e n ). 



(11) 



(12) 



with the dispersion of the incoming electrons' e n . The 
chemical potential /i c will also depend on the scattering 
state as the electrons incoming from the left/right have 
a different chemical potential: ep ± eV/2, respectively. 

Ultimately we are interested in the current through the 
system: 



— ^l(z)-^-ip a (z), 
m oz 



(13) 



FIG. 2. (Color online) The differential conductance G(V) 
normalized to Go = G(V,g = 0), as a function of the bias 
voltage for different widths L of DW in units of the carrier 
wave length Xf- The negative differential conductance occurs 
around £,-. 



for spin a electrons. We are interested in the steady state 
case where there is no charge build up in the wire. The 
current is then 



1 = £&> =- y E-*td*GZ„(z,z';t,t+)\ I , =I . (14) 



We can substitute in our perturbative expression Eq. (|9|) , 
and Eq. (|11[) and work in the limit of zero temperature. 



IV. RESULTS 

The results arc plotted in Figs. [2] and [3j the nega- 
tive differential conductance is marked by a dip in the 
current around the lowest resonance Eq t = e r . We con- 
sider magnetic semiconductor wires^— where the DW 
width can be closer to the Fermi wavelength, the width 
of the wall can be atomically sharp when constrictions 
are presenfe^Zr— . Hence we take \f ~ L, J = 0.005ef- 
The mode energy is e r = 2irJa cs / La, where a cs is the 
cross section and a is the lattice spacing. To keep this 
energy within experimentally reasonable sizes, we con- 
sider quasi-one-dimensional wires where a cs ~ L 2 . This 
gives us e r w 0.03ef^, and for the unrenormalized cou- 
pling we take go = 1QJL At low temperatures in such 
systems electron-correlation effects tend to renormalize 
the coupling strengths! 29 ' 40 Hence we take the effective 
coupling as g e ff ~ 3o£ 1_7 i with £ ~ (3 J = 5 and 7 = 0.8. 
The scattering strengths from the DW will be similarly 
renormalized. 

Of experimental relevance is the differential conduc- 
tance which can be defined as G(V) = dl/ dV. This is 
plotted in Fig. [5] for different DW widths, and the nega- 
tive differential conductance occurs around the resonant 
energy level e r . Higher modes of the DW's motion would 
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FIG. 3. (Color online) The current I(V) as a function of the 
bias voltage. I is measured in units of To = I{eV = £f,<? = 
0). The inset shows the total current, where the feature is 
dimly visible. In the main figure we plot the second order 
correction to I(V): 0(g 2 ). Here the drop in current due to 
the resonant mode coupling is visible. Note however that 
of course the differential conductance is the experimentally 
relevant quantity. 



couple similarly at different energies for larger applied 
bias. 

We note that the model developed in this work is also 
relevant to interacting quantum dot systems involving 
regions with a noncollinear magnetization ! 30 ' This is 
the subject of ongoing work. 

With regard to an experimental realization we note 
the following: NDMR is sizable for DWs in magnctic- 
semiconductor-based nanowires such as those realized in 
Refs. I34H361 . We developed the above theory however 



for electrons as carriers whereas in III-V compounds the 
carriers are holes. As follows from the above treatment, 
the underlying mechanisms for the emergence of NDMR 
are however the DW scattering and inferences of the car- 
riers coupled to the modes of the DW. These elements 
will also be present for carrier holes despite their more 
complicated electronic dispersion, and hence we expect 
NDMR to also be present in this case. 

V. SUMMARY 

In summary, we have calculated the current through a 
ferromagnetic wire with a domain wall present, or equiv- 
alcntly a thin ferromagnetic strip between ferromagnetic 
wires of opposing magnetization. Particular attention 
has been paid to the effects of the lowest modes of the 
DWs motion caused by the presence of this current. In 
order to solve for the density build up around the do- 
main wall subject to an applied potential we used a sol- 
uble model for a sharp DW. The negative differential 
magnctorcsistance expected from the excitation of these 
modes by the current electrons is clearly visible in the 
differential conductance depicted in Fig. [5] This qualita- 
tively new effect offers new opportunities for applications 
in spintronics and domain wall logics^ similar to those 
that utilize NDR in charge-based electronics. 
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